AULA 


Variação de parâmetros 


META: 

Descrever o método conhecido por variação de parâmetros usado 
para resolver equações diferenciais ordinárias lineares não homogêneas. 
Além de apresentar o oscilador harmônico. 

OBJETIVOS: 

Ao fim da aula os alunos deverão ser capazes de: 

Aplicar eficazmente o método de variação de parâmetros. 

Resolver problemas referentes ao oscilador harmônico. 
PRÉ-REQUISITOS 


Os conhecimentos da aula 6. 
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8.1 Introdução 


Caros alunos essa é a nossa oitava aula. Na aula de hoje a- 
presentaremos mais um método de resolução de equações lineares. 
Esse método é destinado a resolver equações lineares não homogênea, 
ele é um pouco mais abragente do que o aprendido na aula anterior 
pois possibilita resolver equações lineares não homogêneas com co- 
eficientes não constantes quando se conhece uma solução geral da 


equação homogênea, associada. 


8.2 Resolvendo equações lineares não homogêneas. 


O método que apresentaremos aqui é conhecido por variação de 
parâmetros. 


Considere a E.D.O. linear 
ant) +: + a(a)y! + ao(x)y = g(x), 


onde an(x),--- ,ao(x) e g(x) são funções contínuas numm intervalo 


I, com an(x) £ O em à. Quando n = 1 obtemos 
y a P(a)y = f(x), 


onde P(x) = ag(x)/a(x) e f(x) = g(x)/a(x). Vimos que a 
E.D.O. homogênea associada y' + P(ax)y = O tem solução geral 
dada por ye(x) — ay(x) — cre! Pldr. Para acharmos yp(x) 
usaremos o método de variação de parâmetros, o qual consiste 
em encontrar uma função uj(x) tal que yp(x) = u(z)yi(x), onde 
y1 é solução da E.D.O. homogênea associada, seja uma solução 


particular de y' + P(x)y = f(a). 
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Substituindo yp(x) = u (x)y (x) na E.D.O. temos 8 


let) + Pla] = H(0) & met ta foh + Plan] = 1(2), 


Como y é solução da E.D.O.homogênea a última igualdade se 


reduz a 


Ássim, 


Logo, a solução geral é 
y(a) = cre JS Plxjde ar e) Pl)da / fHa)el P(a)dr o. 


Observe que essa expressão é a mesma obtida quando resolvemos, 
na aula 4, a E.D.O. linear de primeira ordem utilizando fator in- 


tegrante. 


Quando n = 2, temos 


y" + P(o)y' + Q(x)y = f(x), 


onde P(x) = ar(x)/az(x), Q(x) = ag(x)/az(a) e f(x) = g(x:) /az(a). 
Suponha que yi (x) e y>(x:) formam um conjunto fundamental para 


a E.D.O. homogênea associada 
9" + Pl(a)y' + Q(x)y = 0. 


Assim, procedendo analogamente ao caso anterior, para acharmos 
uma solução particular, y, da E.D.O. não homogênea tentemos 


encontrar funções uj(x), uz(x) tais que yplx) = u(a)y(x) + 
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us(x)ya(x) seja a solução particular procurada da E.D.O. não ho- 


mogênea y” + P(x)y' + Q(x)y = f(x). Derivando yp, obtemos 


, Al 


Yp = Un + uoyo + uy + usy 


Yo = uy + 2uíy + 2usys + uy] + usyo + usyo. 
Sustituindo essas expressões na E.D.O. não homogênea, temos que 


Yo + Pl(x)yp + Q(x)yp = f(x) 


é equivalente a 


uty + Zur + Zuoyo + uy + uzyo + uoyy + P(x)Jurgh + ug + uíyy + usga] 


+Q(x)[uign + uay2] = f(x). 


E, organizando, temos 


uy + P(x)y + Q(u)yn] + uslys + P(x)ys + Q(a)yo] + P(x)[uíy + usgo] 


[ZA Es , Hm 4 , / A 1 (4 
+uty + uiyi + usgyo + usgyo + us + vago — f(x). 


Finalmente, 


P(x)uiy + ubyo] + E lui + ubgo] + uy + ubyo = fa). 
(8.47) 
Observe que para obtermos a última igualdade usamos o fato de 
y1 e y2 serem soluções da E.D.O. homogênea associada, ou seja 
w + Plax)y + Q(x)y =0eyy + P(x)yo + Q(x)yo = 0. 
Como queremos achar duas funções desconhecidas uy e uz pre- 
cisamos de duas equações que as relacionem. Observe que podemos 


obter essas duas equações tomando como hipótese que 


uy + usyo = 0. 
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Dessa maneira, a equação (8.47) se reduz a 8 
us + usyo = f(x). 


Portanto para acharmos as duas funções w e u2 basta resolver o 


sistema 


uy +usyo =0 


uy + uso = Fa). 
O sistema acima pode ser resolvido pela regra de Cramer, a qual 
nos dá 


/ / 


b) 


Wo 
W 


onde W é o Wronskiano das funções y e y2, ou seja 


rest y Ya 
NY 
E, 
(0) (0) 
Wi = det a , Ya = det = 
fax) yo w Fx) 


Integrando as expressões obtidas para u/ e u5, obtemos as funções 


u1, uz procuradas. 


Exemplo 8.1. Resolva a E.D.O. 
y +4y = sec(x), -n/2< 1 < 7/2. 


Lembre que antes de qualquer cálculo a equação deve ser 
colocada na forma padrão, como esta já está na forma 
padrão, prossigamos com os cálculos. 


Sabemos que uma solução geral para essa E.D.O. é da forma 
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y(x) = vela) + yp(x), onde y. é à solução geral da E.D.O. ho- 
mogênea associada e yp uma solução particular da E.D.O. não 
homogênea dada. Começemos com a parte homogênea. 

Uma vez que a E.D.O. homogênea associada y” + 4y = 0 tem 
coeficientes constantes, podemos resolvê-la achando as raízes da 


equação auxiliar 


m2 +4 = 0. 


Vê-se que as raízes dessa equação auxiliar são complexas, as quais 
são m, = 2i,mo = —2%. Assim, a solução geral para a E.D.O. 


homogênea, é 


Vela) = cycos2x + cosen 2x. 


Para resolvermos a parte não homogênea, como a função indepen- 
dente sec x não é uma função polinomial nem nenhum dos tipos 
exigidos pelo método dos coeficientes a determinar, temos que usar, 
neste caso, variação de parâmetros. Dessa maneira, procedamos 
como foi explicado acima. Suponha que a solução particular procu- 
rada seja da forma yp(x) = uy + u2y2, onde y1,y> são as funções 
que aparecem na solução geral da E.D.O. homogênea, a saber 
y = cos2r e y = sen2x. Ássim, derivando yp e substituindo 


na equação, obtemos o seguinte sistema 

u/ cos 2x + ussen 2a = 0 

—2uí sen 2x + 2uí cos 2x = seca. 
Daí, 


Wi —sen2Zrsecr |, W5 cos 2x sec 4 
= Us = = 
, 2 


W 2 W 2 


E 
AS 
I 
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uma vez que 


sen2x 


secar 2cos2r 


cos 2x 0 


—Isen2x secr 


—sen2xsecz, 


= cos2x sec x 


AULA 


= 2cos” 2x + 2sen? 2x = 2. 


0) 
Wi — det 
W5 — det 
e 
cos 2x sen 2x 
W = det 
—Isen2Zr 2cos2x 
Portanto, 
K sen 2x sec a 
ui = 
2 
e 


2 1 
= [ SCE [ 5 cosadr À 


1senZg 


da = [sen vdx = cost 


2 cost 


Portanto, uma solução geral da E.D.O. dada é 


1 
5 2senx-—ln|sec x-+tan z|). 


1 
y(a) = cy cos 2x+casen 2x-+cos 2x cos r+sen 2x[sen x — ain sec x+tan |]. 


8.2.1 Equações de ordem superior 


O método também se aplica a E.D.O. linear de ordem superior 


gO + Pra(o)y D+... + Plo)y + Pola)y = f(x). 


Se vela) = ya) ++: + chynla) é à solução geral para a E.D.O. 


homogênea associada, então yp(x) = uy +: +unYn. As funções 


upi=1, 


;n são obtidas integrando as soluções do sistema 


/ 
uy + 


Pad 
UV 


LASRIA 
Ut Te 


p. (n—1) 


Um +ec+ruyn 


RA 
+unyh 


(n—1) 


= fa). 
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Pela regra de Cramer, temos que 


W:; 
us E W' =lsn 
onde 
ty O) Yn 
y 0 y 
W; = det y A i=1,sn 
; (0) : 
Da = 
gel) f(x) ye 
1 
i- égima coluna 
e, 
mn Yi Yn 
, / / 
y Y; y 
W=det| | k Lista 
= = + 
gr!) ru) e 


Exemplo 8.2. Resolva a E.D.O. 


UMA 


y" + 4y' — sec 2x. 


Uma vez que a E.D.O. está na forma padrão, procedamos como 
no exemplo anterior, calculemos uma solução geral para a E.D.O. 
homogênea associada y” + 4y' = 0. A equação auxiliar nesse caso 
é 
3 o 
mº + 4m = 0, 
cujas raízes são my = 0,ma = 2i,m3 = —2%. Assim, a solução 


geral é dada por 


Vel) = cy + co cos 2x + cgsen 2x. 
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E por variação de parâmetros tentemos obter uma solução particu- 8 
lar yp para essa E.D.O.. Suponha que yp = u1y1 +uoy>-+Husy3, onde 
yi = 1,y92 = cos2x,y3 = sen2x. Neste caso, temos que resolver o 


seguinte sistema 


uí + us cos 2x + ugsen 2x =) 
O — 2ussen2Zr + 2uscos2r | =0 
O — 4us cos 2x — 4ussenZr =sec2a 


Ássim, uma vez que 


0) cos 2x sen 2x 1 (0) sen 2x 
Wi = 0 -92sen2x 2cos2x |, Wa =|0 (0) 2cos2x |; 
sec2r —4cos2r -—4sen2r O sec2x —4sen2gx 
e, 
1 cos 2x (0) 1 cos 2% sen 2x 
Wa=|0 -2sen2x 0 W=|0 -Isenda 2cos2 
O -—4cos2x sec2a O -—4cos2zx —4sen2a 
temos que 
Sá =” “ei a e —1/4,u x fame, 


Integrando, obtemos 


1 1 
ui = girl sec2x + tan 2x], uy = —x/4, ug = gil cos 2x]. 


Assim, uma solução geral da E.D.O. dada é 


1 1 
y(x) = ci+c> cos 2x+casen 21+ gil sec 2x-+tan 22] = cos 2x+ sen (2x)In| cos 24]. 
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8.3 Modelagem matemática em E.D.O. lineares 
de ordem superior com coeficientes cons- 


tantes 


8.3.1 O oscilador harmônico 


Estudamos na aula 1 que a equação diferencial que modela o movi- 
mento de um corpo de massa m presso a uma mola considerando 


a inexistência de forças de amortecimento e de forças externas é 
L=——g4, 
mM 


onde k é a constante de elasticdade da mola e m é a massa do corpo 


presso a mola. Aqui trataremos desse problema com mais detalhes. 
Oscilador harmônico livre não amortecido. 


Consideremos um corpo de massa m presso a uma mola, cuja cons- 
tante de elasticidade é k, livre de amortecimentos e forças externas. 
Se o corpo ao ser presso na mola deslocou-a s unidades de compri- 
mento, segundo a lei de Hooke sabemos que a mola exercerá uma 
força restauradora proporcional ao deslocamento sofrido e oposta 
ao mesmo, ou seja, Frest = ks. Ao deslocar a mola s unidades de 
comprimento o corpo se encontra parado e dizemos que o sistema 
massa-mola se encontra em equilíbrio, ou seja, a força peso é igual 
a força restauradora da mola, mg = ks. 

A força exercida pela mola terá sinal negativo, pois ela age no 
sentido oposto ao movimento. Convencionaremos que os desloca- 
mentos medidos abaixo da posição de equilíbrio serão positivos e 


os acima serão negativo, como mostra a figura abaixo. 
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E 


POSIÇÃO DE EQUILÍBRIO 


s 


Figura 8.1: Sistema massa-mola. 


Num certo instante de tempo t o corpo é puxado para baixo dis- 
tendendo a mola x(t) unidades de comprimento, a qual foi solta 
em seguida. Pela segunda lei de Newton esse movimento obedece 
a seguinte equação mg — k(s + x) = mk, onde mg é a força peso 
e k(s + x) é a força restauradora da mola quando esta é disten- 
dida de s unidades (quando o corpo foi presso à mola, lembre que 
ks = mg) mais x unidades quando o corpo foi puxado para baixo. 
Assim, obtemos a seguinte equação diferencial 


dx 


map =my-k(aL+s)=-kz 


ou podemos escrever a equação que modela o movimento do corpo 


livre de amortecimento e forças externas da seguinte maneira 


dz 
de + wlg = 0, (8.48) 


onde w? = k/m. Esse movimento é chamado movimento har- 
mônico simples ou movimento livre não amortecido. Ob- 
serve que a equação (8.48) é uma equação de segunda ordem do 
tipo linear homogênea. Assim, ao resolvermos a equação auxiliar 
encontramos duas raízes complexas dadas por my = wi,m2 = —wi. 


Portanto, a solução geral de (8.48) é 


a(t) = cy coswt + cosenwt, w = Vk/m. (8.49) 
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Observe que x(t) é uma função periódica, de período T = &r e, 


EEE ne” 
consequentemente, sua frequencia é f = 1/7 — 5. (observe que 


como w2 = k/m deveríamos ter w = +,/k/m matematicamente 
falando, contudo w tem um significado físico e, por isso, adotare- 
mos w sempre positivo.) 

Se ci = 0 ou cz = O podemos facilmente traçar o gráfico da solução 
obtida e, dessa maneira, podemos extrair mais informações sobre a 
solução. Contudo, se c, £ 0 e co + O extrair algumas informações 
se torna um pouco mais difícil. A fim de contornar tal fato, pode- 
mos escrever a solução (8.49) de uma outra forma. Dessa maneira, 


afirmamos que a equação (8.49) pode ser escrita na forma 
a(t) = Asen (wt + wo), (8.50) 


onde A = cj + c5 e senwo = 4, coswo = &. De fato, sabemos 
que Asen (wt + wo) = A[sen wt cos wo + sen wo coswt]. Assim, a 


fim de que 
ci coswt + cosenwt = A[sen wt cos wo + sen wo cos wt] 


temos 


Acoswoy=c2, Asenwç = cc. 


Logo, da relação acima concluímos que as constantes 4 e wo são 


Amplitude de movimento tais que A= ci+ciesenuo= 2 Cos ug = É. 

do corpo presso à mola é a 

maior distância do corpo a Portanto, da expressão (8.50) é fácil ver que a amplitude de movi- 

sua posição de equilíbrio. . , , . 
mento do corpo presso à mola é 4. Concluímos que o movimento 
de um corpo de massa m presso a uma mola de constante elás- 


tica k livre de amortecimento e sem força externa é periódico, 


de período T = 27 ,/m/k e a amplitude máxima desse corpo é 
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A= Vea, ou seja, um corpo nessas condições tenderá a os- 
cilar indefinidamente periodicamente. Abaixo apresentamos o grá- 
fico que descreve o movimento do corpo ao longo do tempo, onde 
a amplitude de movimento é 5 e período igual a 27. O gráfico da 
direita está defasado 2 unidades de tempo, ou seja, wo = 2. No 


gráfico da esquerda wo = O 


O fo ç 4 : 4 [6 10 
| | Ea | 
| | | É | [ ] 
| | | | | 
| | | NH | | | 
| | | | | | 
L4 [ | | | / 
Figura 8.2: Movimento Figura 8.3: Movimento 
livre não amortecido livre não amortecido com 


defasagem de 2 unidades de 


tempo. 


Considerando o gráfico da esquerda como exemplo, observe que 
os pontos onde o gráfico corta o eixo horizontal (eixo do tempo 
t) nos mostram os momentos em que o corpo de massa m passa 
pela posição de equilíbrio, « = 0, os pontos do gráfico que estão 
abaixo do eixo t (ordenada negativa) são os pontos cujas orde- 
nadas nos dão a posição do corpo acima da posição de equilíbrio, 
são as posições correspondentes a compressão da mola e, os pontos 
do gráfico que estão acima do eixo t (ordenada positiva) são os 


pontos cujas as ordenadas nos dão a posição do corpo abaixo da 
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posição de equilíbrio, são as posições correspondentes à distenção 


da mola. 


Oscilador harmônico livre amortecido. 


Nesse caso consideramos a situação anterior sujeita a um amorte- 
cimento seja ele causado pelo vento, pelo atrito com alguma super- 
fície ou por outro motivo qualquer. Na física as forças de amorteci- 
mento que atuam em um corpo são, em geral, proporcionais a uma 
potência da velocidade instantânea. Aqui, consideraremos que a 


força de amortecimento é proporcional a potência unitária da ve- 


locidade, ou seja, a força de amortecimento tem a expressão E, 


onde y é uma constante. Assim, a equação de movimento para 


esse caso é 


mÉ a e dx 
de? — Va 


(O sinal negativo na força de amortecimento ou atrito deve-se ao 
fato que forças de atrito são forças contrárias ao movimento.) As- 
sim, organizando a equação acima obtemos a seguinte equação, a 
qual modela o movimento do corpo com amortecimento mas livre 


de forças externas 


2 
EAR 


Quo = di, 
ET É +wlz=0,20=9/m, wº = k/m. (8.51) 


Esse movimento é denominado por movimento livre amorte- 
cido. Observe que essa E.D.O. é também uma E.D.O. linear 
homogênea e, desta maneira, sua solução geral é obtida por meio 


do cálculo das raízes da equação auxiliar 


mê? +2ôm + w? = 0, 
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a qualnos dám; = —-0+v62 — wlemo = —d—v62 — w2, Assim, e | 


2 


teremos três tipos de solução, a depender se 62 — w2 é maior, igual 


ou menor do que zero. Comecemos com o caso 62 — w? > 0. 
CASO [6º -w?>0 
Seô2-w?>0Otemosm=-04+V62-wlems=-0-v62 — w? 
números reais e, nesse caso, pelo exposto na aula anterior, a solução 


geral é dada por 


a(t) = e Hc e Put 4 coevPuét, 


Esse caso é chamado movimento livre superamortecido. 


Figura 8.4: Movi- Figura 8.5: Movi- Figura 8.6: Movi- 


mento livre super- mento livre super- mento livre super- 
amortecido y(x) = —amortecido y(x) = amortecido y(x) = 
e Tre, 4e* — 368, -p +20, 


CASO 11:60? —-w? = 0 
Se 62 — w2 = 0, temos my; = ma = —6. Assim, pelo conteúdo da 


aula anterior, a solução geral é 
nt) = get + cate *. 


Esse caso é chamado movimento criticamente amortecido. 
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Figura 8.7: Movi- 
mento livre critica- 


mente amortecido 


x 


ylr) =e "+ xe”. 


E, finalmente 


Figura 8.8: Movi- 
mento livre criti- 
camnete | amorte- 
cido  y(x) = 


4e* — 3xe””, 


CASO 11:62 —- w2? < 0 
Se 92 — w? < 0, temos my = -d+ivw2? — d2,mo = d— iv'w2 — 62. 


Assim, pelo conteúdo da aula anterior, a solução geral é 


Figura 8.9: Movi- 
mento livre critica- 
mente amorte- 
cido  y(a) = 


—e" +Igxe 


v(t) — e ey cos Vw? — d2t + cosen Vw? — 82t]. 


Esse caso é chamado movimento subamortecido. 


Figura 8.10: Mov. 
livre subamortecido 
y(a) = eT*[—- cos4x— 


sen 4a]. 


Figura 8.11: 
Mov. livre sub- 
amortecido y(x) = 


e? [cos 4x + 2sen 4x). 
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Oscilador harmônico forçado. 8 


Consideremos agora uma força externa agindo sobre o corpo de 
massa m. Aqui cabe observar que se esta força for a favor do 
movimento ela entrará nas equações com sinal positivo, caso con- 
trário o sinal será negativo. Podemos, por exemplo, pensar nessa 
força externa como sendo uma força periódica que uma pessoa faz 
puxando o corpo para baixo toda vez que ele passar pela posição 
de equilíbrio. De qualquer forma, seja qual for a força externa 
exercida, a equação diferencial de movimento do corpo se con- 


siderarmos além da força externa, a força de amortecimento será 


ms — —ka E + f(t). Organizando a equação, obtemos 
diz da 
S+%— +wlr=F(t 
qo +20 + ta = P(b), 


onde 20 = y/m, w2 = k/m, F(t) = f(t)/m. Se o sistema sofrer 
a ação de uma força externa, mas for livre de forças de amorteci- 
mento, teremos a equação 

d2x 


Note que as equações acima são do tipo linear de segunda ordem 
não homogênea. Assim, se a função F(t) for do tipo polinomial, 
exponencial ou das formas senax, cosax ou combinação desses, 
podemos usar coeficientes a determinar para achar uma solução 
geral, caso contrário a única opção que resta é tentarmos resolver 


por variação de parâmetros. 
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8.4 Conclusão 


Da aula de hoje, concluímos que podemos resolver quaisquer 
equações diferenciais lineares não homogêneas, bastando para isso 
conhecermos a solução geral da equação homogênea associada. 
Para acharmos a solução particular da E.D.O. não homogênea 


usamos um método conhecido por variação de parâmetros. 
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RESUMO 


Na aula de hoje aprendemos o método conhecido por variação de 


parâmetros, cujo objetivo é encontrar uma solução particular da 
E.D.O. linear não homogênea. As equações lineares estão ligadas 
a vários problemas reais, nesta aula também estudamos detal- 
hadamente o oscilador harmônico como uma aplicação para tais 


equações. 


PRÓXIMA AULA 


Até agora só aprendemos a resolver E.D.O.ºs lineares homogêneas 
com coeficientes constantes. Em nossa próxima aula veremos como 
resolver um tipo especial de equações lineares com coeficientes va- 


riáveis: a Equação de Cauchy-Euler. 


ATIVIDADES 


Atividade. 8.1. Resolva cada equação usando variação de parâme- 


tros 


a) y" +y = secr 
by! + y = cos?x 


e2z 


o) y" — 4y = 
d) y' +2y — 8y = 20 —e-* y(0) =1,4(0)=0 


1H 


yu" +y =tga. 


Atividade. 8.2. Uma mola está presa ao teto. Quando uma 


N 


massa de 30kg é atada à mola, esta distende-se em 12 em. A 
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massa é removida e uma pessoa segura a extremidade da mola e 
começa a balançar para cima e para baixo com um período de 1 


segundo. Qual o peso dessa pessoa? 


Atividade. 8.3. Um peso de 0,5 kg é atado a uma mola de 1,5 
m de comprimento. Na posição de equilíbrio, o comprimento da 
mola é de 2,48 m. Se o peso for suspenso e solto a partir do re- 
pouso de um ponto 2 m acima da posição de equilíbrio, encontre o 
deslocamento x(t) se é sabido ainda que o meio ambiente oferece 
resistência numericamente igual à velocidade instantânea. Clas- 
sifique o sistema em amortecido, sub ou super amortecido. Dê a 


amplitude desse movimento. 


Atividade. 8.4. Uma massa de 2 kg é atada a uma mola cuja 
constante vale 32 N/m. Uma força igual a f(t) = 68"? cos 4t atua 
no sistema a partir de t = O. Encontre a equação de movimento 


na ausência de atrito. 


Atividade. 8.5. Uma massa de 0,5 kg é atada a uma mola que 
tem constante de elasticidade igual a 6 N/m. A massa parte do 
repouso 2 m abaixo da posição de equilíbrio e o movimento subse- 
quente está sujeito a uma força de amortecimento igual a metade 
da velocidade instantânea. Encontre a equação de movimento se 


o peso sofre a ação de uma força externa igual a f(t) = 10 cos3t. 
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